A detailed analysis of the spin-crossover reaction of H2S binding to heme and the six-coordinated FeP(Im)-HS- porphyrin complex.

The potential energy surfaces of the H2S binding to iron-porphyrin (FeP) with the imidazole (Im) ligand via intersystem crossings are investigated by using density functional theory. The minimum energy intersystem crossing point (MEISCP) between the quintet and triplet states (MEISCPTQ) for the Fe(II)P(Im)-H2S complex is located at a Fe-S distance of 3.39 Šwith only 1.1 kcal/mol above the quintet state minimum. The second spin-crossover point, where a change from the triplet to the singlet state occurs, comes at a much shorter Fe-S distance of 2.79 Å, and the MEISCPST is located at 3.7 kcal/mol above the triplet state minimum. The nature of the chemical bonding along the Fe-S reaction coordinate from the ground state singlet to the quintet state along the path to the separated species is analyzed. An inspection of the vibrational modes reveals that the largest contribution to the triplet-quintet transition around the quintet and triplet state minimum comes from the symmetric shrinking of the pyrrole units of the porphyrin ring, indicating that the related reaction coordinate plays a main role in the intersystem crossing. The fully optimized structures of the Fe(II)P(Im)-HS- complex corresponding to three different spin multiplicities (M = 1, 3, 5) are characterized by a bent Fe-H-S conformation. The binding of the hydrosulfide anion to Fe(II)P(Im) in the quintet state induces a 0.2 Šdisplacement of the Fe atom out of the nitrogen porphyrin (Npyr) plane. The fully optimized structure of the ground state of Fe(II)P(Im)-HS- agrees well with experimental data for the corresponding heme models.

Modeling the hydrogen sulfide binding to heme.

The binding of hydrogen sulfide to a model heme compound is investigated by coupled-cluster singles-doubles augmented by a perturbative triple excitations, CCSD(T), and density functional theory, DFT. The minimum energy path for the H2S addition to an isolated heme center of the heme protein is evaluated by adopting as a model the heme compound FeP(Im) (P = porphyrin; Im = imidazole). The FeP(Im)-H2S aduct is bound by 13.7 kcal/mol at the CCSD(T) level of theory. Relaxed potential energy curves for the lowest lying spin states of the H2S to FeP(Im) binding using DFT reveal that the binding process is associated with a "double spin-crossover" reaction with the existence of long-distance van der Waals minima only 5-7 kcal/mol above the FeP(Im)-H2S ground state. The fact that the energy of the singlet ground state of FeP(Im)-H2S is so close in energy to the dissociation products FeP(Im) + H2S points towards the reversibility of the H2S adsorption/desorption process in biochemical reactions.

Theoretical study of nitrodibenzofurans: A possible relationship between molecular properties and mutagenic activity.

In this study we present a theoretical investigation of the molecular properties of nitrodibenzofurans (NDFs) and dinitrodibenzofurans (DNDFs) and their relation to mutagenic activity. Equilibrium geometries, relative energies, vertical ionization potentials (IP), vertical electron activities (EA), electronic dipole polarizabilities, and dipole moments of all NDFs and three DNDFs calculated by Density Functional Theory (DFT) methods are reported. The Ziegler/Rauk Energy Decomposition Analysis (EDA) is employed for a direct estimate of the variations of the orbital interaction and steric repulsion terms corresponding to the nitro group and the oxygen of the central ring of NDFs. The results indicate differences among NDF isomers for the cleavage of the related bonds and steric effects in the active site. The results show a good linear relationship between polarizability (<α>), anisotropy of polarizability (Δα), the summation of IR intensities (ΣIIR) and the summation of Raman activities (ΣARaman) over all 3N-6 vibrational modes and experimental mutagenic activities of NDF isomers in Salmonella typhimurium TA98 strain. The polarizability changes with respect to the νsNO+CN vibrational mode are in correlation with the mutagenic activities of NDFs and suggest that intermolecular interactions are favoured along this coordinate.